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ABSTRACT 

The presence of narrow high-temperature emission hnes from nitrogen-rich 
gas close to SN 1987A has been a principal observational constraint on the 
evolutionary status of the supernova's progenitor. A new analysis of the complete 
five-year set of low and high resolution lUE ultraviolet spectra of SN 1987A 
(1987.2-1992.3) provide fluxes for the N VA1240, N IV] A 1486, He IIA1640, 
O III] A1665, N III] A1751, and C III] A1908 lines with significantly reduced 
random and systematic errors and reveals significant short-term fluctuations in 
the light curves. The N V , N IV] , and N III] lines turn on sequentially over 15 to 
20 days and show a progression from high to low ionization potential, implying 
an ionization gradient in the emitting region. The line emission turns on 
suddenly at 83 ± 4 days after the explosion, as defined by N IV] . The N III] line 
reaches peak luminosity at 399 ±15 days. A ring radius of (6.24 ±0.20) x 10^^ cm 
and inclination of 41?0 ± 3?9 is derived from these times, assuming a circular 
ring. The probable role of resonant scattering in the N V light curve introduces 
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systematic errors that leads us to exclude this line from the timing analysis. A 
new nebular analysis yields improved CNO abundance ratios of N/C= 6.1 ± 1.1 
and N/0= 1.7±0.5, confirming the nitrogen enrichment found in our previous 
paper. From the late-time behavior of the light curves we find that the emission 
originates from progressively lower density gas and that the emitting region 
has a multi-component density structure. We estimate the emitting mass near 
maximum (~ 400 days) to be ~ 4.7 x 1O~"^M0, assuming a fiUing factor of unity 
and an electron density of 2.6 x 10^ cm~^. These results are discussed in the 
context of current models for the emission and hydrodynamics of the ring. 

Subject headings: Stars: Circumstellar Material — Supernovae: General — 
Supernovae: Individual(SN 1987A) — Ultraviolet: Spectra — Methods: Data 
Analysis 

1. INTRODUCTION 

Narrow emission lines from SN 1987A were first seen with the International Ultraviolet 
Explorer {lUE ) satellite in 1987 May (Pransson et al. 1989, hereafter Paper I). The 
lines detected were N VAA1239, 1243, N IV] A1486, He IIA1640, O III]AA1661, 1667, 
N III] AA1747-1754, and C III] AA1907, 1909. The absence of the usually strong Si IV A1400 
and C IV A1550 resonance lines can be attributed to absorption of these transitions by 
the same ions in hot interstellar gas in the LMC. The presence of N V in the SN 1987A 
spectra, on the other hand, is in fact an indication that there is a negligible column density 
of interstellar N V in the LMC. A remarkable feature of the early observations of SN 1987A 
was that the N III] , N V , and C III] lines increased in fiux linearly with time. The He II , 
N IV] , and O III] lines were weak and it was difficult to say anything conclusive about their 
light curves. The increase in the fiuxes stopped just after the last observations discussed 
in Paper I, ~ 400 days after explosion. Subsequent observations showed that the C III] , 
N III] , and N V lines reached a maximum ~ 405 days after the explosion (Sonneborn et al. 
1988). After this epoch the observations reported in this paper show a steady decay of the 
lines. Earlier versions of the lUE emission line light curves, but not the new measurements 
discussed in this paper, have been presented by us in several papers (Sonneborn et al. 1990; 
Panagia et al. 1991; Sonneborn 1991) and an independent analysis by Sanz Fernandez de 
Cordoba (1993). Preliminary versions of some of the new light curves have been shown by 
Pransson & Sonneborn (1994) and Plait et al. (1995). In a complementary paper to this one. 
Pun et al. (1995) analyzed the UV spectrophotometric evolution of the SN 1987A debris 
from day 1.6 to 1567 in lUE spectra covering 1150-3300A. They found good agreement 
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with UV continuum observations of SN 1987A with Hubble Space Telescope {HST ) and 
ground-based photometry where they overlap. 

The presence and temporal development of these highly ionized lines have been 
explained as a result of recombination, cooling and light echo effects (Lundqvist & Pransson 
1987; Chevaher 1988; Paper I; Lundqvist & Pransson 1991, henceforth LP91). The initial 
pulse of soft X-rays accompanying the breakout of the shock wave ionizes the circumstellar 
gas expelled by the progenitor star. As the outburst evolves an increasingly large volume 
of the ionized gas becomes visible because of light travel-time effects. If the recombination 
time is long compared to the epoch of observation, the observed flux is expected to grow in 
proportion to the emitting volume within the light echo paraboloid. For a non-recombining 
ion emitting at constant temperature this leads to a flux increasing monotonically with time, 
until the whole structure is within the light echo paraboloid. In reality, the temperature 
and ionization decrease with time, causing the evolution to become highly complicated. 
Detailed modelling of the early line emission showed that the peak radiation temperature at 
the time of the breakout was between 4 x 10^ K and 8 x 10^ K (Pransson & Lundqvist 1989; 
Lundqvist 1991). A subsequent analysis, taking the ring geometry into account, resulted in 
a peak radiation temperature in excess of 10^ K (Lundqvist & Pransson 1996, henceforth 
LP96). These emission line observations therefore provide important information about the 
first moments of the supernova explosion. 

Another important piece of information was the large nitrogen enrichment derived 
from the emission line ratios, primarily the N III] / C III] and N III] / O III] . Abundance 
ratios of N/C = 7.8 ± 4 and N/0 = 1.6 ± 0.8, by number, found in Paper I strongly suggest 
CNO enriched gas. This has proven to be one of the main observational constraints on 
the progenitor evolution (e.g., Saio, Nomoto, & Kato 1988; Woosley 1988; Weiss 1991; 
Podsiadlowski 1992), implying that the progenitor was in a post He-core burning phase at 
the time of the explosion. 

After the discovery of the UV emission lines, circumstellar emission has also been 
observed at other wavelengths. In the optical a large number of lines from both medium 
ionized species and low ionization ions have been detected (Allen, Meikle, & Spyromilio 
1989; Wamplcr k Richichi 1989; Wampler, Richichi, & Baade 1989; Crotts & Kunkel 1991; 
Kahn & Duerbeck 1991; Meikle et al. 1991; Menzies 1991; Wang 1991; Gumming 1994). 
Extended loop-like structures surrounding SN1987A were first imaged in Hq;-|-[N II] and 
[O III] emission by Crotts, Kunkel, & McCarthy (1989). Even in these early images of the 
inner emission region appeared to be a hollow oval. These structures were also seen in 
optical continuum, indicating the presence of time-delayed reflections by dust. Imaging 
with the New Technology Telescope resolved the inner emission region of radius ~ 1" and 
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several outer structures (Wampler et al. 1990; Wang & Wampler 1992). HST images 
in [O III] AA4959, 5007 (Jakobsen et al. 1991, 1994; Plait et al. 1995) revealed the now 
well-known thin ring geometry of the inner region, rather than a spherical shell as had been 
previously hypothesized. Further HST observations showed that the outer structures are 
concentrated mainly to two off-center rings, which extend ~ 2'.'5 to the north and south 
of the supernova (Burrows et al. 1995). Crotts, Kunkel, & Heathcote (1995) show that 
the toroidal geometry extends also to regions of the circumstellar medium exterior to the 
ring, and that there may be some connection between the inner and outer rings. They 
also find that the inner ring is circular to a high degree of accuracy. In addition to the 
geometry, the dynamics and origin of the inner ring is severely constrained by the emission 
line widths. In Paper I we showed that the C III] AA1907, 1909 lines were very narrow 
{FWHM < 30 km s^^) and unshiftcd with respect to the LMC ISM toward the supernova 
(v = 285 ± 6 km s^^). Crotts & Heathcote (1991) and Gumming et al. (1996) found that 
the ring expands with a velocity of 10.3 ± 0.4 km s~^. They, as well as Wood & Faulkner 
(1987) and Mcikle et al. (1991), also find a velocity gradient over the nebula, which is a 
strong indication of a ring geometry rather than a spherical shell. 

Prompt radio emission from SN 1987A, the result of shock interaction with the 
progenitor's stellar wind close to the photosphere, was observed during the first few days 
after explosion (Turtle et al. 1987). This emission decayed away after the first week. 
However, around day 1200 there was a revival of the radio emission and a steady increase 
in the flux was observed (Staveley-Smith et al. 1992, 1993) which continued through 
day 2750 (Ball et al. 1995). Chevaher (1992) and Duffy, Ball, & Kirk (1995) attribute 
the resurgence of the radio emission from SN 1987A to the interaction of the supernova 
ejccta with the density enhancement expected at the position of the termination shock in 
the blue supcrgiant wind of the progenitor. Further evidence of circumstellar interaction 
is the increasing latc-timc X-ray emission observed by ROSAT (Beuermann, Brandt, & 
Pietsch 1994; Gorenstein, Hughes, & Tucker 1994; Hasinger, Aschenbach, & Triimper 1996). 
Chevalier & Dwarkadas (1995) point out that the density of the shocked blue supergiant 
wind is insufficient to explain the late radio and X-ray fluxes, as well as the apparent 
slowing down of the supernova shock (Staveley-Smith et al. 1993). They propose that the 
supernova shock is now moving into the denser {ue ~ 10^ cm~^) shocked red supergiant 
wind which was ionized by the progenitor. 

In this paper we discuss the complete set of UV emission hues in SN 1987A derived 
from lUE low and high dispersion SWP spectra from the explosion up to 1992 April 20 
(day 1882). A new method of analysis has allowed us to reduce systematic and random 
errors in the line measurements compared to Paper I and the light curves of Panagia et al. 
(1991) and Sanz Fernandez de Gordoba (1993). The evolution of the line fluxes, leading to 
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times for turn-on and maximum of the line emission, enable us to determine the absolute 
dimensions of the emitting gas independent of the imaging observations. The implications 
of these results in connection with ground-based and HST observations are discussed. In 
§2 we discuss the observations, in §3 the measurements of the line fluxes, and in §4 the 
main observational results are given. The paper concludes with a discussion (§5) of ring 
geometry, abundances, density, emitting mass, and comparison with existing ring models in 
hght of the observational data. 

2. OBSERVATIONS 

The ultraviolet spectra presented in this paper were obtained with the lUE SWP 
camera between 1987 March 14 and 1992 April 20. All spectra were taken in the 10" x 20" 
short wavelength large aperture (SWLA), the minor axis of which is approximately parallel 
with the spectrograph's dispersion direction. The unique characteristics of the SN 1987A 
field demanded special care during the acquisition process to ensure accurate target 
centering in the SWLA, as discussed below. The observations were carried out at both 
/C/E Observatories, located at NASA/Goddard Space Flight Center and ESA's Villafranca 
ground station (VILSPA). The SWP image number, time of observation, and exposure 
duration for each low dispersion spectrum are listed in Table 1. The observation date refers 
to the starting time of the exposure, expressed in days after outburst. The outburst time is 
assumed to be that of the 1MB neutrino detection (1987 February 23.316 UT, Bionta et al. 
1987). 

Twenty SWP high dispersion spectra of SN 1987A were obtained during the UV 
emission line phase (75 < t < 1800) days. These observations and their analysis are 
discussed below in §4.4. 

Following the rapid decrease of the supernova's UV flux shortward of 2000 A in the 
first week of the outburst, the continua of two early-type stars appeared in SWLA spectra 
of SN 1987A (Kirshner et al. 1987). The identification of the progenitor hinged on close 
scrutiny of these two previously unknown stars in the LMC located within arc seconds of 
the supernova's position. Analysis of their spatial separation in SWP low-dispersion spectra 
taken several weeks after outburst showed that Sk —69° 202 had disappeared (Gilmozzi et 
al. 1987; Sonneborn, Altner, & Kirshner 1987). These same two stars, referred to as star 
2 and star 3, after their position in the tabulation of West et al. (1987), have dominated 
UV fiux in SWP large aperture exposures of SN 1987A since 1987 March. Stars 2 and 3 are 
offset from SN 1987A by 2"91 (p. a. 318°) and 1'.'63 (p. a. 118°), respectively, as measured 
in recent HST WFPC2 images (Romaniello, 1996, private communication. These updated 
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positions are within the measurement errors of Walborn et al. (1987). 

Several known sources contribute to the observed spectrum in the lUE large aperture 
exposures taken at the position of SN 1987A. In addition to the supernova, there is the 
emission from circumstellar material (the subject of this paper), the continua of stars 2 and 
3, and a very faint star located about 9" from SN 1987A which enters the end of the SWLA 
every six months when the aperture has the appropriate orientation. These contributions 
were taken into account in the measurement of emission line strengths. Stars 2 and 3 are 
the most important secondary sources for purposes of this paper. 

The shape of emission lines and supernova continuum are affected by the detailed 
structure of the combined stellar spectra. Correcting the observed SWP spectrum requires 
more than the subtraction of a stellar background spectrum. By virtue of its location 
near the south ecliptic pole, the LMC is observable 365 days per year and the position 
angle of the lUE large aperture changes by about one degree per day. As the orientation 
of the SWLA rotates on the sky during the year, the projection of stars 2 and 3 on to 
the spectrograph dispersion axis also varies. Stars 2 and 3 appear to rotate about the 
supernova's position (relative to the aperture) and their projection onto the direction of 
spectral dispersion has a 365-day periodicity. The spatial separation (4'/5) corresponds 
to ~ 1200 km s~^ (~7 A) when the stars are aligned parallel to the dispersion direction. 
The principal spectral features in the stellar spectra shift in wavelength because star 2 
contributes about 70% of their combined flux below 2000 A. 

Accurate centering of the supernova in the SWLA was necessary in order to obtain 
reliable emission line light curves. During the time interval spanning the observations, 
SN 1987A covered over 13 magnitudes of visual brightness. The lUE Fine Error Sensor 
(FES) was used to acquire the supernova, providing position measurements of its center of 
light from day 1 (V~ 4.5), through visual maximum (day 86, V~ 3.0), and up to day ~ 885 
(V~ 13.4). After that time the supernova was acquired as a blind offset from the nearby 
star Sk -69° 203, 138'/5 almost due N of SN 1987A (p.a.=2?63). The coordinates (B1950.0) 
used for the bhnd offset maneuver are a = 05*^35^" 48!8, 5 = -69°15'39'.'4 for Sk -69° 203 
and a = 05^35"^ 50!0, 5 = -69°ir58'.'0 for SN 1987A. The offset maneuver was tested and 
cahbrated between days 700 and 850 when SN 1987A was still bright enough (12 ^ y ^ 13) 
to directly acquire it with the FES. The measured maneuver errors on these occasions were 
±1 FES unit in each axis, or ±0'.'3 total pointing error on the plane of the sky. We therefore 
believe that there were negligible differences in aperture centering over the five years of the 
observations which might affect the emission line spectrophotometry. 
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3. DATA ANALYSIS 

The procedures used to measure the emission Une strengths in low dispersion spectra 
are described in this Section. All line fluxes were measured relative to the net supernova 
spectrum, that is the observed spectrum with the contribution of stars 2 and 3 subtracted. 
This procedure, which in principle is the same as used in Paper I, assumes that the relative 
contributions of stars 2 and 3 to the observed spectrum are constant. This will be true in 
general provided that the supernova is centered in the SWLA so that both stars are also in 
the aperture and that the stars themselves are not highly variable in the ultraviolet. This 
last point is addressed in §4.1. The principal steps in this process are outlined below. 

3.1. SPATIALLY RESOLVED SPECTRAL DATA 

Our analysis uses the spatially resolved spectra from the extended line-by-line (ELBL) 
data files produced by the standard lUESIPS processing (Turnrose & Thompson 1984, as 
updated by Munoz Piero 1985 for ELBL processing). The ELBL file is a two-dimensional 
spectral image extracted from the lUE photometrically corrected image along lines of 
constant wavelength. For the SWP camera it contains about 700 points in the wavelength 
domain (1150-1980 A, ~ 1.2 A/pixel) and 110 points in the spatial dimension (0'.'68 per 
linear pixel). 

Specialized interactive software was written to correct ELBL pixels corrupted by 
energetic particle hits and bright spots. For each image the area inside the extraction slit 
and in the background areas were inspected for bright spots and cosmic ray effects. When 
found in the background extraction slit, corrupted pixels were replaced with with data 
values corresponding to the average of immediately adjacent pixels. Extreme care was 
exercised when correcting bright spots close to the spectrum. Within the extraction slit 
bright spots were replaced with the average of part of a single row of pixels (parallel to the 
dispersion) at the same distance above or below the dispersion line. 

The spectrum was extracted from the corrected image using a rectangular extraction 
slit, similar to the standard lUESIPS extraction procedure described by Turnrose & 
Thompson (1984). The spatial height of the rectangular spectral extraction slit for the SN 
1987A spectrum was 24 lines (16"). For comparison, the lUESIPS standard slit heights 
are 18 and 30 lines for point source and extended source extractions, respectively. The 
extraction slit height of 24 lines was chosen to include the wings of the point spread function 
of stars 2 and 3 at their maximum spatial separation in the spectrograph format, i.e., when 
their position angle was perpendicular to the dispersion direction. The customized slit 
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height minimizes the background contribution to the extracted spectrum. The standard 
point source extraction sht is too short while the standard extended source extraction sht 
is longer than necessary. Use of a slit-weighted extraction technique, such as 'optimal' 
(Kinney et al. 1991) was not used because of the variable, non-point-source nature of these 
spectra. 

The extracted spectra were corrected for two calibration changes. First, a secular 
shift in the wavelength scale between 1984 and 1988, and second the long-term decline in 
SWP sensitivity. Starting in 1984 a systematic error developed in the SWP low dispersion 
wavelength scale and grew linearly with time. The error resulted from a slow divergence 
of true dispersion relation from the extrapolation of the mean time-dependent dispersion 
constants. New dispersion constants were implemented in lUESIPS in 1988 April at 
GSFC and 1988 September at Vilspa (Thompson 1988). The SWP low dispersion spectra 
processed between 1987 February and 1988 September contained an error in the wavelength 
calibration amounting to ~ 2.5 A. Thompson (1988) provided an algorithm to fit the 
time-dependent error in the wavelength scale. This relation was used to correct all SWP 
low dispersion spectra of SN 1987A effected by this calibration change. 

Long-term changes in the sensitivity of the lUE instrument were compensated for by 
using the correction scheme of Garhart (1992a). Between 1225 and 2000 A the SWP low 
dispersion sensitivity has decreased linearly since 1979.5 at wavelength-dependent rates of 
—0.2% to —1.2% per year. The Garhart sensitivity correction method gives results similar 
to Bohlin & Grillmair (1986) for the period in which they overlap. The extracted fluxes 
were also corrected for the temperature dependence of the SWP camera response (—0.5% 
per °G, Garhart 1992b). Finally, the spectral fluxes were placed on the lUE absolute flux 
scale using the 1980 absolute cahbration (Bohhn & Holm 1980). 

3.2. REFERENCE SPECTRUM OF STARS 2 AND 3 

There was a period of a few weeks in 1987 when the SWP large aperture spectrum 
of SN 1987A was due almost entirely to stars 2 and 3. This occurred between mid-March 
and early April, before the emission lines appeared. The SWP large aperture spectra taken 
during this period were examined to determine which ones could be combined to form a 
reference spectrum of stars 2 and 3 . Criteria included no overexposure, good signal-to-noise 
at wavelengths shortward of 1500 A (i.e., texp ^ IGOmin), and no apparent signal from 
SN 1987A, particularly longward of 1800 A. Three spectra were selected: SWP 30592 
(day 27.3, t^^p = 300 min), SWP 30637 (day 33.1, t^^p = 185 min), and SWP 30743 (day 
44.2, texp = 225 min). These spectra were reduced in the same manner described in §3.1. 
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SWP 30637 and 30743 were registered in wavelength with respect to SWP 30592 by cross 
correlation. The mean spectrum is the average of the three weighted by the individual 
exposure times. The resulting SWP spectrum of stars 2 and 3 is shown in Figure 1. This 
spectrum is similar, though not identical, to that used by Pun et al. (1995) for their study 
the UV continuum properties of the SN 1987A outburst. They derived their star 2 and 
3 UV spectrum by averaging late-time lUE spectra (1300 < t < 1500 days), when the 
emission lines were extremely weak, in order to recover the UV spectrum of the supernova 
early in the outburst. In contrast, we use early-time spectra for stars 2 and 3, before the 
emission lines appear, so that we can follow the decline of the emission lines at late times. 
Our spectrum of stars 2 and 3 is in good agreement with the FOS spectrum of star 2 (Fig 1 
of Scuderi et al. 1996), given the differences in spectral resolution and calibration between 
the two instruments. Star 2 contributes ~70% of the flux in this wavelength region. 

3.3. BACKGROUND STAR SUBTRACTION 

Subtraction of the reference spectrum of stars 2 and 3 from the observed spectrum 
requires that the two wavelength scales be registered to correct for any offset. Such a 
wavelength shift could be due to either an offset in the center of light from the center of the 
SWLA or a shift in the spectral features of stars 2 and 3 in the observed spectrum. In the 
case of the former, a I'.'O pointing error in the dispersion direction corresponds to a spectral 
shift of ~ 300 km s~^, or ~ 1.6 A. The wavelength shift required to properly align stars 
2 and 3 with the observed spectrum was computed with an interactive cross-correlation 
program from the NASA/GSFC lUE Data Analysis Center software hbrary. 

Figure 2 shows two SWP spectra with aperture position angles differing by 156°. 
The shifts in the observed spectra, relative to the reference spectrum of stars 2 and 3, 
are evident, as is the difference in the widths in the strong stellar hnes. Relative to the 
reference spectrum, the spectra are shifted by —459 and +243 km s~^, as determined by 
the cross-correlation analysis. Figure 3 shows the wavelength registration shifts for the 
SN 1987ASWP data set (Table 1). There is a ~ 1 year periodicity in the data, indicating 
that the position of stars 2 and 3 relative to the spectrograph-aperture geometry is a 
major source of this variation. Within about ±0'.'5 (~ ±150 km s^^) the wavelength 
registration shifts are dominated by target acquisition and pointing errors. The net 
spectrum is computed by subtracting stars 2 and 3 from each observed spectrum after 
shifting the spectrum of stars 2 and 3 by the amount determined by the cross correlation 
and interpolating it to the wavelength scale of the observed spectrum. 
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4. EMISSION LINE RESULTS 

4.1. FLUXES AND LIGHT CURVES 

The flux in each emission hne was measured in the net spectrum where local continuum 
determination is simplified for most of the emission lines of interest. The subtraction of 
stars 2 and 3 removes strong stellar and interstellar absorption features, including Si II 

A1260, Si III A1300, C IIA1335, and Al III A1670. The wavelengths defining the center 
and red and blue edges of each line in each spectrum were determined interactively. After 
subtracting stars 2 and 3 the local continuum was assumed to be fiat, except where the 
SN UV fiux is strong. The latter is confined to A 1700 A and t ^ 700 days. The line 
strength was integrated by a trapezoid algorithm. An example is shown in Figure 4 where 
the observed and net spectra near the N V line in SWP 36801 (day 894) are shown on the 
same flux scale. The N V hne is easily distinguished in the net spectrum. 

Table 1 contains the fluxes for the N VA1240, N IV] A1486, He IIA1640, O III] A1665, 
N III] A1750, and C III] A1908 lines measured in 119 spectra taken on 106 dates. The 
1225 - 1980 A history of SN 1987A between days 85 and 1820 is shown in Figure 5 where 
the observed net fluxes have been averaged into nine time intervals. The development of 
the six emission lines is clearly seen. In addition to these emission lines, other features are 
evident in the net spectra. The supernova continuum longward of 1700 A is prominent 
until day ~ 700 (see also Pun et al. 1995). There are additional features which rise and 
fall in the same time period between 1280 and 1540 A. Of particular interest is the narrow 
feature at 1400 A, which we tentatively identify with O IV] AA 1397- 1407. Although we have 
no identifications for the other features in the spectra, they are believed to be part of the 
Fe II blanketed continuum of the supernova. The formation of such "pseudo-emission" 
features in the early spectra of SN 1987A is discussed by Cassatella et al. (1987) and Lucy 
(1987) and at late times by Fransson (1994) and Li & McCray (1996). Local maxima occur 
where line blocking is less important and the photons can escape, i.e., at wavelengths where 
there are gaps in the list of strong lines. Similar pseudo-emission features occur in classical 
novae. Hauschildt et al. (1994) modelled and discussed the formation and evolution of 
these features in Nova Cygni 1992. 

The resulting light curves for the N V , N IV] , He II , O III] , N III] , and C III] lines 
are shown in Figure 6. The light curves show the observed line fiuxes, without correction 
for interstellar extinction (see §5.2). It was sometimes difficult to determine the continuum 
level at the base of the C III] emission line, even after subtracting stars 2 and 3, because 
of the UV fiux from the supernova debris increased in this region of the spectrum (see 
Fig. 5). A check of the accuracy of the low dispersion C III] flux measurements was made 
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by comparing them with C III] fluxes from high resolution observations where continuum 
determination relative to the narrow emission is much easier. The high dispersion data 
for SN 1987A are discussed below (§4.4) and the corresponding C III] fluxes also shown 
in Figure 6d. The good agreement between these independent measurements, well within 
the scatter of the low resolution data, illustrates the overall accuracy of the continuum 
determination. The new emission line light curves are clearly superior to those published 
previously (Paper I, Sonneborn et al. 1990, Panagia et al. 1991, and Sanz Fernandez de 
Cordoba (1993) and may effect results derived from them, for example Gould (1994, 1995). 

The size of the lUE large aperture prompts us to consider the degree to which emission 
from the outer rings might contribute to the UV light curves, particularly at early times 
(say, t < 600 days). If one adopts the geometry suggested by Burrows et al. (1995) the 
rising branch of the light curve of the southern outer ring (SOR) should extend from 
about one month after explosion to about three years later, whereas the northern outer 
ring (NOR) should have become visible about two years after outburst. In this geometry 
the NOR would not have been entirely visible in the earliest images (Crotts, Kunkel, & 
McCarthy 1989) but may have been "complete" by late 1989 in the ESO images (Wampler 
et al. 1990). The UV light curves do show a weak secondary maximum around 1000 
days, consistent with the epoch when both outer rings would have reached their maximum 
emission (Panagia et al. 1996, in preparation). Therefore, emission from the NOR would 
not contribute to the early lUE observations, and emission from the SOR would have 
contributed less than 1/3 of its total strength. Present day HST imaging shows that the 
NOR, the SOR and the inner ring have comparable total emission in optical lines. However, 
the average density in the inner ring is higher than that in the outer rings so the time decay 
of the inner ring emission would be faster than that of the outer rings. As a consequence, 
the total emission of the inner ring at the time of its maximum must have been considerably 
higher than that of the outer rings combined. In fact, ESO observations around the time of 
the tentative secondary maximum show that the outer rings contribute less than 10% of the 
total [O III] emission (J. Wampler & L. Wang 1992, private communication). We therefore 
believe that the contribution of the outer rings to the lUE emission hue hght curves is 
neghgible. 

Figure 4 shows that the subtraction of stars 2 and 3 is not perfect. Due to fixed-pattern 
and other small-scale noise, the scatter in the net spectrum near the emission line limits the 
accuracy of continuum subtraction, and hence the significance of short-term fluctuations in 
the emission line light curves. The accuracy of the continuum subtraction is quantified by 
the parameter 5Fc, the local continuum residual (Eqn. 1), which is the difference between 
the net spectrum and the adopted local continuum, averaged over a bandpass comparable 
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to the width of the emission hne, 

SF, = {FnetiX) - FiocaliX)) d\. (1) 

5Fc was measured near the N V hne because the adjacent continua of stars 2 and 3 have 
strong stellar and interstellar features, which if not subtracted accurately will add significant 
and erroneous structure to the net spectrum. Perfect continuum subtraction would yield 
5Fc = 0. 6Fc was computed in a bandpass centered at 1262 A with a full width of 10 A, 
the characteristic width of the N V emission line at zero intensity. Figure 7 shows the N V 
light curve (same data as Fig. 6a) and SF^. The error bar icr = 1.82 x 10~^^ erg cm~^ s~^ 
corresponds to the measurement uncertainty of individual data points for the six emission 
lines, indicating that the continuum subtraction procedure works well over the entire 5-year 
span of the data set. 

The accuracy of the continuum subtraction is insensitive to the known variability of 
the two stars. First, a change in the continuum flux results in a vertical shift of the net 
spectrum. The continuum in the net spectrum is constant (i.e. horizontal in Fig. 4), except 
in the presence of the strong SN continuum, so that the emission line is measured relative to 
this shifted continuum level. Second, it is now well estabhshed that only star 3, a Be star, 
shows any significant variability. Walborn et al. (1993) showed that star 3 has photometric 
variability ranging from ±0.15 mag in U to ±0.3 mag in R and /. Since star 3 contributes 
only ~ 30% of the UV flux, this level of variability is within the measured value of a. 

The measurement uncertainty of the line fluxes implies that the signiflcant short-term 
fluctuations in the light curves are astrophysical in origin, and are not instrumental effects. 
The increased scatter in the N V and N III] light curves between day 900 and 1100, in the 
N IV] line before day 400, and the large values of N III] and N V near maximum are actual 
fluctuations in the emission line fluxes. We flnd no evidence that these spectra are affected 
by cosmic ray hits or other spurious effects and conclude that the fluctuations are real. The 
largest N V flux point (day 485) remains suspect, however, as it occurs signiflcantly after 
maximum brightness and there is no related response in the other lines. 

4.2. EMISSION LINE MAXIMUM 

The times when the lines flrst appear and when they reach maximum brightness are 
two important light curve parameters used to estimate the size and inclination of the ring. 
With a circular ring geometry the emission lines are expected to appear suddenly as the 
hght echo paraboloid makes flrst contact with the ring (Dwek & Felten 1992). In general, 
the highest ionization stage available should be used to measure these parameters. 
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It is evident from Figure 6 that the six UV hnes reach maximum at approximately 
the same time. However, only the N V and N III] lines have large fluxes and well-defined 
maxima necessary for a more quantitative analysis. The actual times of maxima are 
t^iJJ = 399 ± 15 and t^^^ = 429 ± 19 days (418 ± 12 days if the high data point at day 485 
is excluded), as determined by high-order polynomial fits to the light curves. Comparison 
of the N V and N III] in Figure 8 shows that while their post-maximum decline is very 
similar, the light curves before maximum have significant differences. Figure 8 shows a 
large deficiency of N V relative to N III] between 150 and 400 days and that N V peaks 
about three weeks after N III] . 

The N VA1240 emission from the far side of the ring may be delayed by 1-2 weeks by 
resonance scattering in the ring and the bipolar nebula, as discussed by LF96. However, to 
obtain an even longer delay, and thus later t^^, like the one observed, additional resonance 
scattering is needed. The most hkely source is in the H II region interior to ring discussed 
by Chevalier & Dwarkadas (1995). Preliminary calculations by us show that there may be 
a sufficient optical depth in N V A1240 through this region along the line of sight to the far 
side of the ring. The large deficiency of N V relative to N III] is therefore probably due to 
resonant scattering of N V photons by N"*"^ ions in and outside the ring. This biases the 
N V flux maximum and complicates its interpretation by introducing effects other than 
light travel time and local recombination. In view of the uncertainties in interpreting the 
N V light curve, we believe t'^^x^ provides the best measure of the time of maximum, tmax, 
and we use this in our subsequent analysis. 



4.3. EMISSION LINE TURN-ON 

Figure 9 shows the N III] , N IV] , and N V fluxes up to day 120. All three lines appear 
very abruptly, as predicted. There is also an ionization dependence in the turn-on times. 
N V appeared first, after day 70.2 and before day 80.6. N IV] appeared next, between 
80.6 and 85.2 days. N III] appeared last, between 85.2 and 90.9 days. While it is possible 
that N III] was present at a very weak level between day 70 and 85, the measurements are 
consistent with zero flux based on the scatter before day 60. Since the semi-forbidden lines 
are not effected by scattering, N IV] , the next highest ionization state after N V , should 
closely track the inside edge of the ring. However, as discussed below (see also Fig. 6 in 
LF96), the N V light curve is signiflcantly affected by resonance scattering which introduce 
large systematic errors. Consequently we believe N V is also inappropriate for measuring 
trise- We take the sharp rise in the N IV] flux as the turn-on time and find Uise = 83 ± 4 
days (estimated error). The uncertainty in the turn-on time is set by the spacing of the 
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observations, and consequently trise ma-y not be resolved in the lUE data. 

On the near side of the ring resonance scattering of N V by ions in the ring and 
any other intervening highly ionized circumstellar gas should delay the turn on of this line 
compared to that of other nitrogen lines. Instead, the observations show an earlier turn 
on of N V . The most likely explanation to this apparent inconsistency is the ionization 
stratification of the ionized gas. The effect of stratification on the turn on of N V should be 
larger if there is an H II region interior to the ring from which there also may be weak N V 
emission. Another possible source of early N V emission is the bipolar nebula close to the 
ring, but away from its equatorial plane where the nebula may not be ionization bounded. 
This part of the nebula may have a lower density than the ring causing it to glow early in 
N V , but not in other UV hues due to long recombination time scales. 

4.4. HIGH DISPERSION SPECTRA 

The first SWP high dispersion spectrum of SN 1987A during the narrow emission line 
phase was taken on 1987 Nov. 25 (day 275) and used in Paper 1. Nineteen additional high 
dispersion SWP spectra were obtained between then and day 1652. The observations are 
summarized in Table 2. Using standard lUESIPS high dispersion processing, extracted 
spectra were calibrated using the flux cahbration of Cassatella et al. (1994), which provides 
an accuracy of about 4% with respect to low resolution spectra longward of ~ 1400 A and 
increasing by no more than a factor of two at shorter wavelengths. All the lines seen in 
low resolution were also detected in high dispersion spectra taken near peak brightness. 
However, only the N III] and C III] lines were detected over a sufficiently long time base 
and with sufficient signal-to-noise to be able to discuss their evolution and compare them 
with the low dispersion data. 

Table 2 also contains the fluxes for C III] AA1906.68, 1908.73 and N III] AA1749.68, 
1752.16. The other members of the N III] multiplet were either too weak or, in the case of 
A1748.61, badly effected by a camera rcscau. Figures 10 and 11 show the spectra for these 
lines. The line widths are consistent with instrumental broadening (FWHM ~ 0.18 A at 
N III] and ~ 0.20 A at C III] ) for the time period where hne widths are reliably measured 
(up to day ~1000 for C III] ). 

The relative strengths of the C III] and N III] multiplet components arc well known 
diagnostics of electron density (Nussbaumer & Schild 1979; Nussbaumer & Storey 1979; 
Berrington 1985). In Paper I we showed that rie = 1 — 3 x 10^ cm~^ at the time of the 
first high dispersion spectrum in Table 2 (day 275). The larger number of spectra now 
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available makes it possible to examine the line ratios for evidence of density variations. 
Figure 12 shows the C III] A1907/A1909 and N III] A1750/A1752 line ratios before day 
1000. Calculated /igorZ-^igog intensity ratio as a function of electron density are shown in 
Figure 13 for three temperatures using a 6-level C III atom and recent atomic data. The 
C III] line ratio of 0.8 ± 0.2 up to day 408, implies rie = 4.7 ± 2.0 x 10"^ cm~^, assuming 
Tg ~ 50, 000 K. After ~ 600 days, Iigor/hdod ~ 1-5 ± 0.4. The average electron density 
for the C III] emitting gas decreases by a factor of ~ 3 or more from ~ 5 x 10^ cm~^ 
before ~ 400 days to ^ 1 — 2 x 10^ cm~^ thereafter, depending on Tg. Near ~ 600 days 
Te ~ 30, 000 K may be more appropriate, favoring the lower estimate of Ug- The increasing 
uncertainty in the later data points makes the line ratios after ~ 800 days unreliable for 
density analysis. 



5. DISCUSSION 

The absolute dimensions of the emitting region can be derived from the shape of the 
UV light curves with some basic assumptions about its geometry, mainly a circular thin 
ring. For this purpose we summarize the main features of the recombination/cooling - light 
echo scenario. Hydrodynamic models of the breakout of the supernova shock wave through 
the photosphere of the progenitor show that a strong pulse of soft X-rays is emitted with a 
duration of a few hours (Shigeyama, Nomoto, & Hashimoto 1988; Woosley 1988; Shigeyama 
& Nomoto 1990; Blinnikov & Nadyozhin 1991; Ensman & Burrows 1992). The photospheric 
temperature, and also the radiation temperature, decrease from ~ 10^ K to ~ 2 x lO'' K 
in about 12 hours. The ionizing flux from the supernova is therefore exclusively emitted 
during the first day, and in contrast to the interstellar visual light echoes, we can consider 
the emission from the supernova as instantaneous. The soft X-rays ionize the circumstellar 
gas on a time scale of hours, leading to highly ionized species like the observed N V , or 
even higher (LF91; LF96). At the same time the radiation heats the gas to temperatures 
exceeding ~ 10^ K. The subsequent evolution of the gas, occurring on a time scale of weeks 
to years, is determined by the coupled cooling and recombination of the plasma. 



5.1. GEOMETRY 

The recombination and cooling times are set by the temperature and density, 
tree ~ 1/ (ttrec'^e) , SO gas at low dcusity recombiucs and cools on a longer time scale. The 
finite velocity of light implies that only a region bounded by the light echo paraboloid, given 
by r = ct/(l — cos 6*), can be observed from the Earth at time t. Here r is the distance from 
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the supernova, t the time of the observation, and 9 the angle between the hne of sight and 
the direction to the emitting gas as seen from the supernova. In the case of a spherical shell, 
with radius Rs, the whole structure is seen after a time t„tax = '^Rg/c, while for a circular 
ring tmax = (1 + smi)R/c, where R is the ring's radius and i is the inclination of the plane 
of the ring with respect to the line of sight {i = 0° is face-on). For ions which recombine 
and cool slowly, the observed flux increases nearly linearly up to tmax- If the density is high 
enough recombination occurs rapidly and the light curves of different ions reach maximum 
simultaneously followed by decay on the individual recombination time scales. In the slow 
recombination case the light curves will remain at a roughly constant level after maximum. 

LF91 showed that for a spherical shell the gas was marginally optically thin in the 
continuum. In the case of a ring geometry, however, the gas must have significant optical 
depth in the EUV continuum because of the larger column density needed to produce the 
same emission (LF96). If the ionization of the gas were radially constant (as is probably 
not the case), the light curves of medium ionization ions formed by recombination, like 
N III , would be delayed with respect to those formed by photoionization at the time of the 
outburst, like N V. However, from the outset the ring will be optically thick to ionizing 
photons, resulting in a radial ionization gradient with several ionization zones (Lundqvist 
1992). Medium ionization ions like N III and N IV are created with high abundances 
outside the N V - N VI zone. Models using a ring geometry are discussed in detail by 
LF96, and are consistent with the observed light curves. The delayed turn-on of N IV] and 
N III] (Fig. 8) relative to N V supports this model. 

The HST images clearly resolve the ring-like shape of the emitting gas in line and 
continuum radiation (Jakobsen et al. 1991, 1994; Plait et al. 1995; Burrows et al. 1995). 
While the extent of the region perpendicular to the ring plane remains to be determined by 
modelling of the fines, there are additional indications from the lUE observations for a thin 
ring geometry. The most direct evidence for this comes from the abrupt rise of the emission 
line fluxes. As discussed in §4.1, we adopt the N IV] turn-on time {trise = 83 =b 4 days) as 
the initial rise time. In §4.1 we also determined t^ax = 399 ± 15 days from the N III] line. 
The delay from an inclined circular ring is given by trise = (1 ~ sini)R/c. We therefore have 

R = c{trnax+trise)/2 = (6.24 ± 0.20) X 10^^ Cm, and Sini = c{tmax -trise) /2R = 0.656 ±0.039, 

or i = 41?0 ± 3?9. The inclination has also been measured from the HST FOC images, 
assuming a circular ring. Jakobsen et al. (1991) found i — 43° ± 3° and Plait et al. (1995) 
found i — 44.0° ±1.0°. The agreement between these values shows that the ring is nearly 
circular (see also Gould 1994; Lundqvist 1994; Crotts, Kunkel, & Heathcote, 1995). 

Our value of R has direct bearing on the estimated distance to SN 1987A using the 
method devised by Panagia et al. (1991), and refined by Gould (1994, 1995). Gould 
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finds trise — 75.0 ± 2.6 days and t^ax — 390 ± 1.8 days, which combined with the angular 
size of the ring in [O III] at t > 1200 days from Plait et al. (1995) gives a distance, 
dsN < 46.77 ± 0.76 kpc. Our slightly larger values of trise and t^ax and more realistic errors 
translate directly into a larger value and uncertainty for the distance, dsN = 48.6 ± 2.2 kpc. 
A more complete statistical treatment of the light curves and the time of maximum and 
related systematic errors effecting the distance to the supernova is beyond the scope of 
this paper but is discussed by Lundqvist et al. (1996, in preparation). We emphasize that 
systematic errors may be significantly larger than statistical errors for several reasons. 
For example, the method assumes perfect circularity of the ring. Although this may 
be a reasonable approximation for the overall geometry of the ring, there may be local 
variations of the radius around the ring. In particular, a small inward kink or near-absence 
of gas on the far side of the ring would give a smaller t^ax than for a circular ring, thus 
underestimating dsN- Second, the structure of the [O III] ring at t > 1200 days was 
probably different from the structure of the gas emitting the UV lines at t <^ 400 days. 
Indeed, LF96 found that the ionization structure of the ring changes dramatically in only a 
few hundred days. The UV lines are somewhat more sensitive to temperature and somewhat 
more peaked towards the inner edge of the ring than forbidden line emission (e.g., [O III] ). 
The analyses of Panagia et al. and Gould may be effected by a systematic error arising 
from the assumption that the optical forbidden lines and UV lines come from the same gas 
(but not observed at the same epoch). Depending on the radial distribution of the density 
in the ring, the [O III] structure could either have a smaller, or larger, mean radius than 
the early UV line emission. We also point out that their analyses used UV line fluxes from 
two independent reductions (labeled GSFC and Vilspa shown by Panagia et al. (1991) of 
the same set of spectra, thus the number of independent data points is too large by a factor 
of two. 



5.2. ABUNDANCES 

The increased number of observations compared to Paper I, reduced measurement 
and systematic errors, and improved signal-to-noise ratios make it important to repeat 
the abundance analysis. In order to convert the ionic ratios obtained from the line ratios 
to elemental abundances, assumptions about these relations have to be made. More 
specifically, in Paper I it was assumed that N/C = N III /C III and N/0 = (N III + 
N IV) /O III. The calculations by LF91 and LF96 show that, because of a balance of 
recombinations and collisional ionizations, this is a good approximation until the decay of 
the lines begin. After this epoch the gas temperature is too low for collisional ionization 
to be important for C III, N III, and O III, and the lines decay on their respective 
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recombination time scales. These are quite different for the ions involved, and consequently 
the relations above break down. Before ~ 200 days the lines were weak and the fluxes 
therefore uncertain. Consequently, in this analysis we restrict ourselves to 200 < t < 500 
days. The collision strengths are taken from Berrington (1985), Keenan et. al. (1986), 
Blum & Pradhan (1992). 

The extinction toward SN 1987A has been discussed in several papers and there is now 
general agreement. Walker & Suntzeff (1990) have measured the reddening in the field 
of SN 1987A by studying the {U — B) and {B — V) colors of 23 early-type stars in the 
vicinity of the supernova. They determine the reddening to be E{B — V)totai = 0.17 ± 0.02. 
Fitzpatrick & Walborn (1990) obtained optical and UV spectra of the LMC BO. 7 la 
supergiant ~2' from SN 1987A (cf. §2). They conclude that E{B - V)totai = 0.16 for 
SN 1987A, with Galactic and LMC components of E{B - V)Gai = 0.06 {Ax from Savage 
& Mathis, 1979) and E{B - V)lmc = 0.10 {Ax for 30 Dor from Fitzpatrick, 1985). We 
adopt the Fitzpatrick & Walborn extinction toward Sk -69° 203 for SN 1987A. The total 
extinction correction factors for the six emission lines are listed in Table 3. In Paper I we 
adopted E{B — V)totai = 0.20, the extinction correction for which is also given in Table 3 
{E{B — V)Gai = 0.06 and E{B — V)lmc = 0.14). The line luminosities are computed by 
correcting for extinction and the distance to the supernova. Assuming isotropic radiation 
and a distance of 50kpc, the dilution factor is 3.0 x 10^^. 

The ratios depend weakly on the temperature, and in Table 4 we give the resulting 
values of N/C and N/0 as function of this parameter. In Paper 1 we argued that 
Te ~ 5 X 10"^ K from [O III] observations (Wampler & Richichi 1989) and the models in 
LF91. Using this temperature we find that N/C = 6.1 ± 1.1 and N/0 = 1.7 ± 0.5, where 
the errors are purely statistical. Accepting an (unrealistically large) uncertainty in Tg from 
10^ K to 10^ K, corresponds to a range in N/C of 4. 7 — 15.7. The statistical errors in the 
N/C ratio are therefore considerably less than the systematic. Because the lower limit is 
N/C = 4.7, and a more realistic value is N/C = 6.1, the conclusions in Paper I on the large 
nitrogen enhancement are unaltered. For comparison, LF96 derive N/C = 5.0 ± 2.0 and 
N/0 =1.1±0.4 from their photoionization/recombination model. The two methods are 
in good agreement. 

There is a weak tendency for the derived N III /C III ratio to decrease with time. For 
400 < t < 800 days we find N III /C III = 5.3 ± 1.5 for = 5 x 10^ K. The decreasing 
N III/C III ratio is consistent with the expectations from the models in LF91, and justifies 
limiting the abundance analysis to t ^ 500 days, when N/C = N III/C III is a reasonable 
approximation. 

The nebular analysis has a number of limitations. A nebular model gives a reliable 
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estimate of the ionization structure under steady-state conditions or immediately after the 
gas has been ionized. As the ionizing radiation vanishes, the ionization state evolves away 
from this structure due to different recombination time scales for the different elements (C, 
N, and O). For example, at tmax the part on the very far side of the ring can be roughly 
explained by a nebular structure, but the gas on the very near side has had tmax — trise 
days to evolve away from such a structure. Because the observed emission at tmax is an 
overlay of emission from different parts of the ring which have had different times to evolve 
off the steady-state structure, a pure nebular analysis is bound to have some systematic 
error. In contrast, LF96 did not have to assume an ionization structure for their models 
since it and the temperature are calculated in a self-consistent manner. This results in 
smaller systematic errors related to ionization equilibrium. Given these uncertainties, the 
CNO abundance ratios in this paper and in LF96 are in good agreement. In the near future 
spatially resolved spectroscopy of the inner ring with STIS on HST will make possible this 
type of nebular analysis as a function of position around the ring. 



5.3. EMITTING MASS 

iFiom. the He IIA1640 line we can estimate the mass of the He II emitting gas, i.e., the 
mass of He III zone. In order to avoid having to compensate for light travel time effects, we 
evaluate the He II flux near day 400 so that the full volume is observed. Near tmax ±100 
days the observed He II flux was (9 ± 2) x 10~"'^^erg cm~^ s~^ (estimated uncertainty), and 
the luminosity 1.0 ± 0.2 x 10^^ erg s~^, assuming E(i? — V") = 0.16. The He II emissivity 
does not vary much with time (LF91 and Fig. 6e ), so LhcII ~ / jeffiTg) fig nneiii dV , 
where je// is the effective emissivity of the A1640 line. In the He III zone riHeiii ~ nne- 
Furthermore, Wang (1991) finds that UHe/iT-H ~ 0.2, and we find for the total mass of this 
zone ^ 

M(He III) ^ 1.4 X 10-2 ( -) ( '^^) Mq. (2) 

^ ' V2.6 X 10^ cm-3y \5xWkJ ° ^' 

This mass estimate only refers to the He III zone. As models by LF96 show, the total mass 
of ionized gas can be several times larger. A better estimate of the emitting mass is obtained 
from the H/3 luminosity. Wampler & Richichi (1989) estimate -L(H/3) 1.4 x 10^^ erg s~^ 
on day 310 {E{B -V) ^ 0.16). Multiplying this by a factor 3/2, since only ~ 2/3 of the 
total ring was within the light echo paraboloid at that time, and again using nHe/'f^H ~ 0.2, 
we obtain a total mass of the H II region 

M(H H) « 4.7 X 10-2 ( ^ -\ ( '^^] ' Mq (3) 

^ ^ V2.6 X 104 cm-3y V5 x WKJ ° ^ ^ 
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on day 395. Because of the high temperature, both the He II hne and H/3 are affected by 
colUsional excitation. The mass estimates are therefore upper hmits to the actual masses. 
(Sec LF96 for an estimate which takes these model dependent effects into account.) Our 
estimates, however, suggest that the He III zone occupies ~ 30% of the total ionized region. 

The first if^T observations of SN 1987A were made on day 1278 (Jakobsen et al. 1991), 

and we can combine a mass estimate using B.f3 at this time with the observed thickness of 
the ring to infer the angular extent of the ring, A^, perpendicular to the ring plane, as 
seen from the supernova. From the [O II] and [N II] lines the density and temperature of 
the gas were at this time rig ~ 1 x 10^ cm~^ and Tg ~ 1.7 x 10^ K (Wang 1991; Gumming 
et al. 1996). The flux in R/3 on day 1280 was ~ 7.2 x lO^^^erg cm^^ s^^ (Wang 1991), 
and the luminosity 2.2 x 10^^ erg s~^ (for E{B — V) = 0.16). The total mass was therefore 
4.3 x 10-2(ne/10^ cm-=^)-i(Te/1.7 x 10^ K)^ '^^ M©. Although characteristic of a different 
density component, this mass is not very different from that near day 400. Jakobsen et 
al. (1991) estimate the thickness of the ring in the plane of the sky, AR±, to ~ 10% of 
the radius, or ~ 6 x 10^^ cm. Plait et al. (1995) found a slightly larger value, ~ 15%. 
Because later images obtained shortly after day ~ 2500 using the COSTAR corrective 
optics (Jakobsen et al. 1994) tend to favor the lower estimate of ~ 10%, we will use this in 
our analysis. We note, however, that the emitting region could have been somewhat more 
extended at day 1280 without contradicting the observations by Jakobsen et al. (1994). 
Assuming a volume filling factor within the ring, e, we get 

A. « 1.4».-' ( ^) " ( ^' , )°" . (4) 

VI X 10^ cm-3; \^ 0.1 J V1.7x104k; 

This angular extent corresponds to a thickness along the line of sight, it!||. 



0-012e-M— --^ , (5) 



0.91 



R VlxlO^cm-V \^ 0.1 J Vl.7xl04K 

around day 1278, i.e., smaller than the radial extent (~ 0.1 i?), unless e <C 1. For dynamical 
reasons we would if anything expect the radial extent to be smaller than the lateral. A 
solution to this problem may be that the filling factor is considerably smaller than unity, 
and/or that AR±/R ^0.1. The HST [O III] images indeed show considerable structure 
which may indicate e < 1 in the ring. However, it is difficult to see that it can be much 
less than unity if the ring is resolved in the most recent images (Jakobsen et al. 1994), 
unless the ring consists of blobs smaller than the HST resolution, ^O!,'! (~ 7.5 x 10^^ cm at 
a distance of 50 kpc). These conclusions agree with those of Plait et al. (1995) and LF96. 
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5.4. DENSITY COMPONENTS 

The line emission becomes dominated by gas of decreasing density because cooling and 
recombination time scales are inversely proportional to electron density. Therefore, the flux 
seen at early times is likely to come from gas of considerably higher density than at later 
times (LF96), consistent with the order of magnitude decrease in Ue derived in §4.4 from 
high-dispersion lUE spectra. This agrees with optical results for [O II] at day 1280 by 
Wang (1991) who finds Ue ~ 0.7 x 10^ cm~^. A steady decrease in electron density was also 
found by Plait et al. (1995) when they analyzed HST [O III] images spanning the epoch 
1278 - 2431 days. The emitting mass may therefore vary with time. 

Additional information about the density distribution can be inferred from the low 
resolution lUE observations. LF96 find it necessary to include gas with electron densities 

down to ~ 6 X 10^ cm~^ to explain the slow decline of the lines after maximum. Such a 
distribution of densities in the emitting gas is expected from hydrodynamic simulations 
of the circumstellar structure of SN 1987A (Blondin & Lundqvist 1993), which show that 
the highest densities are found close to the ring plane. They also find that there is a large 
amount of low density gas both in and out of the plane, possibly accounting for the gas 
causing the resonance scattering of N VA1240. If the low density gas is ionized up to N V 
and higher, and has a density of ;^ a few x 10^ cm~^, this may also give rise to a weak, 
steady fine emission, as predicted by LF91. They argued that such emission should be 
the sign of the inner parts of the unshocked red supergiant wind. In the interacting winds 
scenario (e.g., Blondin & Lundqvist 1993; Martin & Arnett 1995), the emission could also 
come from the shocked red supergiant wind (LF96) or perhaps the H II region interior to 
the ring (Chevalier & Dwarkadas 1995). The increasing importance of density components 
with He ^ 10^ cm~^ after ~ 1000 days is consistent with the observed decrease in electron 
density inferred from density-sensitive line ratios (C III] , §4.4) and the optical [O II] and 
[S II] lines (Wampler et al. 1989; Wang 1991; LF96). 

There are also some theoretical indications for gas with ng 4 x 10^ cm~^. LF96 model 
three components with different densities and find that the bulk of the C III] emission at 
t ^ 400 days originates in their highest density component (ug ~ 4 x 10^ cm~^), i.e., close to 
what we derive in §4.4. It is likely that there is a continuous range of densities rather than 
components with a few discrete densities. The component with ~ 4 x 10^ cm~^ in LF96 
therefore represents a mean value at the high end of this range, indicating the existence of 
gas with even higher densities. The scatter in the N IV] light curve is considerably larger 
during the first 400 days than at later times (Fig. 6c), which may be a further sign of 
components with rig ^ A x lO"' cm^'^. Because recombination from N IV to N III is rapid, 
and coUisional ionization unimportant, the emissivity of the N IV] line declines rapidly with 
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time at a given point in the emitting gas (see Fig. 2 of LF91). Therefore, the flux of the 
N IV] hnc is not severely smoothed by the light echo, and responds rapidly to changes in 
the density of material being swept up by the light echo paraboloid. When higher density 
gas is encountered, the result is a momentary increase in the observed flux, and may thus 
explain the fluctuations. Fluctuations due to density variations around the ring decrease in 
importance when the full ring is seen and only recombination is important. The narrow 
spikes in N V and N III] between days 350 and 500 (see Fig. 6) are also probably due 
to higher density gas and the 'arcsine' shape of the hght curve, characteristic of a ring 
geometry (e.g., Felten & Dwek 1991; Lundqvist 1991). 

There is also evidence for density structure in the gas after maximum, especially in 
N III] and N V near day 900 to 1000. These fluctuations indicate that there may be regions 
of high density ~ 2 — 3 times further away from the supernova than the ring radius. (This 
flux increase might also arise from the outer rings, as discussed in §4.1.) Several groups 
have detected blobs in the optical and near-infrared which may be responsible for these 
fluctuations. Allen et al. (1989) detected a concentration of He 11.083 fiui emission, and 
they estimate a gas density of ~ 10^ cm~^. Extended He 11.083 fim emission was also 
observed by Ehas et al. (1993). Crotts, Kunkel, & McCarthy (1989) and Hanuschik (1991) 
detected Ha blobs on days 750 and 950, respectively, and Gumming & Meikle (1993) found 
evidence for a high density blob with short life time on day 1344. Gumming & Meikle 
estimate that the gas density in the blob has to be (1 — 2) x 10^ cm~^ and argue that it 
may be accounted for by the high density regions forming off the equatorial plane in the 
hydrodynamic models by Blondin & Lundqvist (1993). 

5.5. COMPARISON WITH MODELS 

Comparison of the observations with a reahstic, self-consistent 
photoionization/recombination model of the circumstellar ring is necessary to test 
the interpretation of these data. Such a model includes solving the coupled, time dependent 
equations for temperature and ionization; separating these two aspects leads to an 
unphysical model. Knowledge of the emissivity as a function of time and position makes it 
straightforward to incorporate the light echo effect. This was first done by Lundqvist & 
Fransson (1987), and later by Fransson & Lundqvist (1989) and LF91. In those models it 
was assumed that the emission originated in a spherical shell of density 2.6 x 10'^ cm~^. For 
the ionizing burst at breakout two SN 1987A explosion models were used, the lOL model 
by Woosley (1988) and the 11E1Y6 model by Shigeyama et al. (1988). Qualitatively, 
both outburst models reproduced the observed relative line strengths surprisingly well. 
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considering the lack of free parameters. The calculations agree fairly well with the rapid 
increase of the N III] and N V lines, and the nearly constant flux in the N IV] line diiring 
the first ~ 400 days. The steady He II flux was found to be a result of collisional excitation. 
The model of LF91 based on Shigeyama et al.'s 11E1Y6 spectrum provided the best 
agreement with both the behavior of the UV lines and the gas temperature, as determined 
from the [O III] AA4363/(4959 + 5007) ratio. The harder burst spectrum by Shigeyama et 
al. was therefore favored. 

A closer examination, however, reveals several deficiencies. Calculations using both 
the lOL and 11E1Y6 models produced too rapid an increase in emission flux, especially the 
N V line, which rises almost linearly with time (Fig. 6a). As discussed by Lundqvist (1994) 
and LF96, this may be the result of resonance scattering of the N V line, not included in 
the early models. LF96 also find that the slow dechne of the N V line may be a result of 
resonance scattering in gas external to the ring. The scattering gas could extend as much 
as ~ 10^^ cm away from the ring plane. The fact that the maxima of the N III] , N IV] , 
N V and C III] lines arc all observed to occur at nearly the same time poses a more serious 
problem for the LF91 models, which in contrast show the N III] and C III] lines peaking 
considerably later than N V . In addition, the optical [O III] lines show a persistent emission 
up to ~ 1500 days, in contrast to the early models. 

There are two important factors not taken into account in the first models which 
probably account for these discrepancies. As acknowledged by LF91, the assumption of only 
one density component is artificial, and the observations clearly indicate that this is not 
true. A combination of density components, or more realistically a continuum of densities, 
ranging from rig ~ 6 x 10^ cm~^ to at least ng ~ 4 x 10^ cm~^, gives considerably better 
agreement, as demonstrated by LF96. However, such a distribution cannot be specified 
in a unique way as it depends on the detailed structure of the ring. The other important 
assumption in LF91 was that of a spherical geometry. A toroidal geometry changes the 
light curve, when integrating over the light echo paraboloid, as discussed by Lundqvist 

(1991) , Luo (1991) and Dwek & Felten (1992). A generic feature of a ring model compared 
to a spherical shell is the rapid rise in the observed fluxes close to maximum light. This 
basically reflects the amount of emitting gas swept up by the light echo paraboloid per unit 
time. Even more important, the larger radial thickness required by a toroidal geometry 
also implies that the ionized gas will be optically thick in the EUV continuum. As already 
explained, the ionization will then be stratified radially, with several distinct ionization 
zones (Lundqvist 1992; LF96). The consequence for the fight curve is to give nearly 
simultaneous maxima for the different ions, as is observed. In addition, Ensman & Burrows 

(1992) find that the ionizing spectra used in LF91 are probably too soft, as was also pointed 
out in LF91. Ensman & Burrows find that the SN spectrum at breakout will have a high 
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energy tail due to electron scattering. Ensman & Burrows also note that light travel time 
effects smooth the time variation of the spectrum, as may also be the case if the progenitor 
was non-spherical (see Lundqvist 1992 for a more detailed discussion). However, LF96 find 
that this effect is unimportant for the ionization of the ring. 

In summary, the light curves are fully consistent with a circular ring geometry, both 
from the rapid rise close to maximum light seen in the N III] and N V lines, the turn-on 
time of these lines, and the near-coincidence of the light curve peaks. A ring geometry also 
makes it easier to understand the low expansion velocity, ~ 10 km s^-*^ (Luo & McCray 
1991; Blondin & Lundqvist 1993). The origin of the ring, however, remains obscure (e.g., 
Podsiadlowski 1992 for a review). The high N/C and N/0 ratios put strong constraints on 
the ring implying that it must have been formed after the dredge up phase. An interesting 
proposal by McCray & Lin (1994) explains the ring as the result of the gradual evaporation 
of a circumstellar disk around the progenitor. They explain the nitrogen enrichment as a 
result of gradual spraying of the inner boundary of the disk by the red supergiant wind. Any 
model for the formation of the ring must also be able to explain the outer rings observed 
by Wampler et al. (1990), Wang & Wampler (1992), Burrows et al. (1995) and Plait et 
al. (1995). The interacting winds model (e.g., Luo & McCray 1991; Blondin & Lundqvist 
1993) may do so, as demonstrated qualitatively by Martin & Arnett (1995), and so may the 
model by Lloyd, O'Brien, & Kahn (1995). The interacting winds model can also account for 
the continuum observations of Crotts et al. (1995) (Martin & Arnett 1995), and the hot gas 
responsible for the resonance scattering of the N V line (LF96). However, the high density 
of the outer rings, a few 10^ cm""^ (estimated by LF96) suggests that the interacting winds 
model needs some adjustment. An interesting extension is that of Chevalier & Dwarkadas 
(1995). They assume a lower mass loss rate during the blue supergiant stage, implying that 
photoionization by the progenitor is important during the shaping of the nebula. A test of 
the models for the formation of the observed structure around SN 1987A will be possible as 
the supernova ejecta strikes the ring shortly after the year ~ 2000 (e.g., Masai & Nomoto 
1994; Luo, McCray, & Slavin 1994; Chevaher & Dwarkadas 1995; Borkowski, Blondin & 
McCray 1996). 
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Table 1. Emission Line Fluxes for SN 1987A Circumstellar Ring 



SWP^ 


Date'' 


t ^ 


V 

corr 


Integrated Flux 


(units: 


1 X 10~^^ 


ere; cm 

o 


-2 s-i) 




(days) 


(min) 


(km/s) 


N V 1240 


N IV] I486 Hen 1640 


O III] 1666 


N III] 1750 


C III] 1908 


30440 


11.3 


300 


-257 


5 


13 


2 




1 




30472 


14.0 


180 


46 


14 


-7 


14 




4 


1 


30512 


18.3 


240 


-46 


2 


-9 


-1 


17 




6 


30522 


19.1 


181 


591 


11 


4 


-2 


7 


16 


-4 


30547 


21.3 


240 


111 


6 




-2 


17 


10 


5 


30592 


27.3 


300 


-125 


3 


7 


19 


-2 


6 


-5 


30637 


33.1 


185 


-106 


2 


10 


18 


-5 


4 


5 


30743 


44.2 


225 


53 


4 


21 


4 


6 


1 


2 


30907 


70.2 


230 


164 


7 


4 


26 


19 


14 


7 


30974 


80.6 


120 


303 


55 


2 


49 




12 


11 


31000 


85.2 


180 


-37 


69 


45 


70 


12 


17 


8 


31040 


90.9 


195 


233 


64 


114 


33 


35 


29 


17 


31064 


96.3 


180 


199 


84 


101 


74 


12 


76 


6 


31125 


105.8 


200 


-7 


105 


77 


55 


27 


117 


13 


31132 


107.8 


195 


42 


92 


95 


67 


42 


109 


6 


31154 


110.8 


180 


-154 


144 


107 


54 


68 


113 


27 


31166 


113.0 


180 


-83 


107 


108 


70 


50 


88 


29 


31245 


121.7 


260 


8 


98 


155 


43 


66 


90 


91 


31273 


128.0 


240 


111 


98 


131 


75 


66 


191 


38 


31319 


138.5 


240 


7 


115 


163 


92 


75 


170 


70 


31371 


147.0 


240 


130 


159 


117 


94 


64 


215 


78 


31420 


154.5 


240 


280 


163 


141 


116 


94 


207 


72 


31462 


160.9 


240 


122 


140 


139 


85 


139 


190 


53 


31534 


168.5 


220 


396 


150 


117 


78 


88 


265 


88 


31592 


177.7 


240 


243 


186 


136 


93 


84 


183 




31651 


185.4 


210 


-282 


146 


92 


95 


141 


222 


94 


31676 


188.8 


240 


-119 


147 


157 


60 


130 


233 


X 


31819/i5 


198.9 


2A0/ 90 


-32/173 


163 


186 


67 


108 


302 


106 


31893 


209.6 


200 


-104 


175 


181 


75 


93 


268 


110 


31954 


220.3 


200 


-926 


131 


148 


74 


99 


265 


101 


32030/5i 


227.7 


240/ 


-U2/-34 


139 


177 


110 


155 


389 


101 


32168 


245.5 


155 


520 










343 


102 
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Table 1 — Continued 



SWP'^ Date"^ texpo^ Vcorr^ Integrated Flux (units: 1 x 10 erg cm ^ s ^) 





(days) 


(min) 


(km/s) 


N V 1240 


N IV] 1486 


Hen 1640 


Oiii] 1666 


Niii] 1750 


C III] 1908 


32219/^0 


252.6 


240/56* 


-352/ -570 


176 


126 


90 


85 


323 


113 


32314 


263.4 


235 


46 


215 


161 


93 


94 


X 


X 


32395 


275.8 


90 


-500 










321 


135 


32404 


276.6 


240 


-98 


234 


221 


96 


97 


X 


X 


32532 


296.2 


240 


-212 


223 


214 


85 


145 


X 


X 


32619/^0 


305.5 


240/50 


-493/ -^55 


233 


202 


77 


132 


344 


177 


32717 


324.2 


190 


-414 


291 


141 




111 


311 


114 


32797 


338.4 


210 


-189 


322 


170 


108 


116 


411 


X 


32879 


349.5 


75 


-45 










426 


162 


32911 


355.1 


212 


-170 


338 


198 


76 


95 


414 


137 


32938 


360.0 


180 


-591 


349 


174 


69 


97 


425 


114 


33035 


375.3 


240 


85 


379 


197 


116 


127 


X 


X 


33105/0^ 


388.0 


240/ 7^? 


b/-211 


375 


211 


114 


126 


482 


206 


33175/75 


400.3 


21Q/90 


1-292 


464 


195 


84 


149 


442 


199 


33280 


416.0 


200 


621 


498 


151 


110 


172 


511 


197 


33331/5^ 


423.2 


240/75 


135/-i55 


383 


202 


80 


173 


444 


181 


33423 


433.6 


80 


-40 










338 


178 


33492 


440.6 


80 


345 










531 


176 


?,U97/96 


441.9 


240/70 


?,/292 


385 


152 


79 


137 


404 


202 


33519/^0 


445.1 


240/5C 


-ns/162 


418 


170 


98 


131 


362 


215 


33644 


458.3 


80 


73 










322 


211 


33725 


471.4 


80 


-326 










324 


143 


33742/^ i 


473.8 


2?,h/80 


lQZ/44 


348 


121 


73 


119 


367 


241 


33799/ 


485.0 


240/50 


U5/153 


561 


133 


85 


114 


327 


178 


33869 


498.7 


195 


-132 


343 


130 


84 


147 


330 


145 


33937 


510.9 


240 


220 


303 


113 


91 


106 


321 


X 


33967/55 


518.7 


251/00 


172/-236 


272 


109 


112 


136 


361 


144 


34057/55 


532.9 


240/50 


172/-183 


314 


118 


95 


88 


266 


127 


34088 


538.5 


200 


347 


267 


141 


75 


99 


294 


120 


34232 


567.5 


200 


115 


245 


105 


86 


92 


270 


103 


34441 


593.5 


150 


-200 


194 








199 


100 


34640 


615.8 


240 


-283 


203 


82 


40 


98 


223 


111 
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Table 1 — Continued 



SWP^ 


Date^ 

(days) 


t ^ 

^GXpO 

(min) 


Vcorr*^ 
(km/s) 


Integrated Flux (units: 1 x 10 


erg cm ^ 

O 


/ 


N V 1240 


N IV] 1486 


Hen 1640 


Oiii] 1666 


Niii] 1750 


C III] 1908 


34751 


632.7 


240 


-400 


197 


61 


85 


72 


222 


100 


34871 


650.7 


163 


-502 


203 






104 


213 


89 


35126 


670.7 


180 


-691 


190 


63 


77 


114 


167 


94 


35308 


686.5 


220 


-332 


238 


74 


65 


67 


184 


99 


35505 


714.3 


220 


-453 


-t CUT 

187 


62 


54 


87 


170 


68 


35686 


740.3 


240 


114 


176 


24 


Ct A 

24 


73 


94 


70 


35822 


755.1 


183 


33 


182 


48 


33 


74 


119 


106 


35940 


771.2 


260 


-29 


152 


55 


55 


50 


112 


66 


36171 


796.2 


210 


59 


135 


63 


40 




123 


74 


36258 


809.9 


222 


136 


159 


21 


45 


87 


131 


49 


36279 


813.2 


260 


267 


143 


68 


54 


55 


129 


50 


36539 


844.0 


260 


96 


134 


25 


48 


52 


99 


45 


36578 


853.8 


310 





142 


35 


51 




106 


54 


36676 


870.9 


260 


141 


151 


28 


47 


87 


101 


33 


36801 


893.9 


260 


93 


129 


18 


44 


55 


120 


58 


36968 


928.5 


270 


89 


137 


7 


28 


44 


67 


43 


37062 


936.5 


230 


-182 


119 


1 


34 


70 


102 


48 


37088 


938.8 


300 


218 


115 


1 


42 


54 


121 


54 


37236 


951.5 


263 


873 


89 


3 


36 


55 


85 


26 


37424 


970.7 


230 


-359 


103 


5 


46 


39 


108 


47 


37574 


992.6 


300 


-501 


122 


9 


51 


55 


91 


52 


37798 


1021.6 


300 


-422 


82 


42 


44 


51 


64 


62 


37973 


1049.5 


285 


-314 


99 


33 


69 


56 


63 


39 


38055 


1061.1 


280 


-748 


92 


-9 


19 


60 


95 


32 


38172 


1081.3 


285 


-332 


101 


4 


11 


42 


69 


78 


38307 


1105.3 


285 


-424 


135 


20 


5 


43 


48 


19 


38336 


1112.0 


275 


-414 


115 


-3 


24 


44 


50 


33 


38536 


1137.2 


285 


90 


89 


7 


4 


57 


51 


34 


38866 


1186.1 


285 


-152 


80 


18 


13 


22 


24 


34 


39300 


1247.0 


285 


-6 


50 


10 


35 


55 


28 


23 


39757 


1317.7 


265 


-287 


90 


-11 


6 


40 


35 


19 


40002 


1344.7 


270 


-550 


65 


-10 


53 


58 


40 


10 
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Table 1 — Continued 



SWP^ 


Date^ 

(days) 


t ^ 

(min) 


V 

* corr 
(km/s) 


Integrated Flux (units: 


1 X 10-^^ 


erg cm ^ 




N V 1240 


N IV] 1486 


Hen 1640 


Oiii] 1666 


Niii] 1750 


Cm] 1908 


40275 


1380.5 


280 


-491 


43 


6 


23 


28 


28 


9 


40858 


1449.3 


280 


-422 


4 


2 


3 


28 


19 




41179 


1490.4 


310 


-110 


17 


-5 


2 


25 


16 


11 


41802 


1566.0 


290 


198 


20 


-2 




20 


10 


3 


42174 


1623.0 


270 


-372 


46 


10 


32 


26 


2 


15 


43049 


1719.6 


340 


-307 


37 


-10 


3 


22 


2 


5 


43393 


1755.6 


340 


-255 


21 


1 


9 


10 


11 


-4 


44017 


1820.4 


340 


-678 


43 


10 


-3 


28 


8 


7 


44142 


1840.3 


360 


-405 


17 


3 


-3 


24 


5 


-1 


44442 


1882.2 


270 


-577 


18 


19 


5 


28 


9 


1 



^SWP image number. On thirteen dates two images of different exposure times were taken 
consecutively to obtain good-quality spectra of all six lines. Data entries corresponding to 
exposures shorter then 100 minutes are given in italics. 

^Observation date, measured in days from the 1MB neutrino event, 1987 Feb. 23.316 (Bionta 
et al. 1987). 

■^Exposure duration in minutes. 
^Cross-correlation velocity correction (see §3.3). 

denotes that the raw spectral data is saturated at this wavelength. 
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Table 2. High Dispersion C III] and N III] Emission Line Fluxes 



SWP 


Date^ 


t ^ 


F 1907^^ 


-^1909'^ 


-^1749'^ 


-f'1752^ 


32394 


275.544 


450 


6.277 


7.085 


15.84 


12.63 


32983 


367.110 


540 


7.513 


9.318 


19.64 


12.72 


33215 


408.167 


880 


7.924 


9.832 


23.76 


15.12 


33536 


448.026 


900 


11.973 


7.996 


22.17 


14.53 


33810 


487.901 


850 


7.139 


5.311 


18.96 


9.84 


34024 


528.874 


952 


6.790 


5.560 


16.25 


12.52 


34420 


590.622 


760 


6.229 


4.578 


14.07 


8.43 


34870 


650.340 


612 


8.094 


5.099 


10.57 


5.72 


35577 


726.232 


925 


5.634 


3.530 


12.69 


6.15 


35929 


771.117 


630 


4.303 


3.468 




6.90 


36399 


834.004 


880 


6.821 


3.311 


9.96 


7.59 


36891 


917.794 


945 


3.234 


2.897 


5.48 


3.10 


37297 


960.701 


975 


3.236 


1.567 


8.05 


3.44 


37805 


1024.533 


970 


2.512 


1.481 


17.34 


3.57 


38230 


1096.065 


610 


2.576 


1.372 


5.42 


2.69 


39093 


1208.039 


965 


1.623 


0.907 


6.66 


5.04 


39644 


1300.799 


930 


1.803 


0.954 


5.98 


3.29 



Table 2 — Continued 



SWP Date^ texpo" -^1907^ -^1909'^ -^1749^^ -^1752^ 

40268 1380.515 1010 1.589 1.071 ■•■ 2.81 

41316 1505.059 625 0.865 0.546 3.27 2.35 

42370 1651.787 960 0.508 0.253 3.81 1.73 

^Exposure mid-point measured in days since 1987 Feb. 54.316 UT 

^Exposure time duration in minutes 

"^Integrated line flux in units of 1 x 10~^^ erg cm~^ s~^ 
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Table 3. Extinction Corrections 



Line E(B-V)=0.16^ Lgs^ E(B-V)=0.20^ Lgg^ 



N VA1240 


5.75 


1.73 


9.32 


2.80 


N IV] A 1486 


4.18 


1.25 


6.19 


1.86 


He IIA1640 


3.76 


1.13 


5.34 


1.60 


III] A1666 


3.68 


1.10 


5.20 


1.56 


N III] A1750 


3.61 


1.08 


5.08 


1.52 


C III] A1908 


3.59 


1.08 


5.01 


1.50 



^Extinction correction for observed flux 

^Linc luminosity in 10^^ erg s~^, per unit flux of 1 x 
2Q-13 gj,g (,j^-2 g-i^ assuming d — 50 kpc and the indicated extinction. 



Table 4. CNO Abundance Ratios 



Te/IO^K N/C Error N/0 Error 



1.0 


15.7 


2.9 


1.4 


0.3 


2.0 


9.4 


1.7 


1.5 


0.4 


5.0 


6.1 


1.1 


1.7 


0.5 


10.0 


1.7 


0.9 


1.8 


0.5 
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Fig. 1. — The reference SWP spectrum of stars 2 and 3 is the exposure- weighted sum of 
SWP 30592, 30637, and 30743. The principal stellar and interstellar absorption features are 
identified below the spectrum. The position of the six UV emission lines are marked above 
the spectrum. A prominent camera reseaux is indicated by 'R'. 
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Fig. 2. — Two SWP spectra of SN 1987A (thicker lines) illustrating the effect of the changing 
orientation of stars 2 and 3 with respect to the dispersion axis. The SWLA position angles 
for the two spectra differ by f56°. The reference spectrum of stars 2 and 3 (cf. Fig. f) is 
shown by the thinner line. 
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Fig. 3. — The wavelength (velocity) registration shift for the SN 1987A SWP database, as 
determined by cross-correlation with the star 2 and 3 reference spectrum, shown as a function 
of time. The right-hand axes show the corresponding spatial and wavelength scales. The 
scatter at the level of ±125 km s~^ is a measure of the accuracy with which the supernova 
is centered in the aperture. The one-year cycle of aperture rotation on the sky produces the 
long-term variation. 
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Fig. 4. — Continuum subtraction for N Vin SWP 36801. The upper solidlme is the observed 
SN 1987A spectrum. The upper dashed line is the reference spectrum of stars 2 and 3 from 
1987 March (cf. Fig. 1). The difference spectrum (^A^ - stars) is shown by the solid line at 
the bottom of the figure on the same fiux scale. The solid horizontal line is the continuum 
level established for this N V fiux measurement [shaded area). The dashed horizontal line 
marks the zero fiux level. The N V emission feature is more easily distinguished in the 
difference spectrum. The dotted lines above and below the difference spectrum are the iblcr 
residuals for the reference spectrum. Cross-correlation techniques were used to register the 
SN 1987A and star 2 and 3 spectra; a +93 km s~^ correction was applied to the reference 
spectrum in this example. 
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Fig. 5. — The mean SWP spectrum of SN 1987A (minus stars 2 and 3) is shown for nine 
time intervals spanning the first 1820 days of the outburst. Each mean spectrum is the 
average of five exposures over the indicated time intervaL They are plotted on the same 
fiux scale and are offset for clarity. The dashed horizontal lines mark the zero level for each 
spectrum. The six emission lines measured are marked at the top of the figure, as is the 
possible identification of IV] multiplet UV 0.01. 
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Fig. 6. — The light curves spanning day 11 to 1882 for the six UV emission hnes in SN 1987A. 
The open circles in Fig. 6d are the combined C III] A1906.68 and A1908.73 fluxes from the 
lUE high dispersion spectra (Table 2). 
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Fig. 7. — The N VA1240 light curve [filled circles) and continuum residual, 5Fc (Eqn. 
1), near the line are shown as a function of time. SFc is a measure of the accuracy and 
consistency of the continuum placement used to compute the flux in the line (cf. Fig. 
4). The results shown for N V are representative of the six lines measured. The error bar 
(iblcr = 1.82 X 10"^'* erg cm~^ s~^) is the measurement uncertainty for individual data points. 
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Fig. 8. — Comparison of the N III] and N V light curves near maximum. The large deficiency 
of N V relative to N III] before maximum is believed to be the result of resonant scattering 
by N"*""* ions in the highly ionized gas close to the ring along the line of sight to the supernova. 
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Fig. 9. — Comparison of the three nitrogen emission hne hght curves before day 120. The 
sudden onset of the hues also indicates a progression from higher to lower ionization potential. 
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Fig. 10. — Line profiles of C III] multiplet UV 0.01 detected in lUE high dispersion spectra 
of SN 1987A between days 275 and 1652. The two vertical dashed lines mark the positions 
of the A1906.68 and A1908.73 components at a heliocentric velocity of +285 km s~^. 
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Fig. 11. — Line profiles of N III] multiplet UV 0.01 detected in lUE high dispersion spectra 
between days 275 and 1652. The five vertical dashed lines mark the positions of the multiplet 
components at a heliocentric velocity of +285 km s~^. 
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Fig. 12.— Evolution of the N III] A1749/A1752 and C III] A1907/A1909 line ratios between 
days 275 and 1000. 
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Fig. 13. — Calculated C III] line ratios (/1907/A909) as a function of Ue for a 6-level C III atom 
are shown for electron temperatures of 30,000 K, 50,000 K, and 80,000 K. The atomic data 
are from several sources: A-values for C III from Nussbaumer & Storey (1978), AUard et al. 
(1990), and Fleming et al. (1994), energy levels are those of Mendoza (1983) and Aggarwal 
et al. (1986), and collision strengths from Seaton (1987), and Keenan et al. (1992). 



